clc; %清楚前面的指令
close all; %关闭matlab所有窗口
clear all; %清空工作区
I=10; %初始患者人数
N=10000; %总人数
S=N-I; %剩余的就是易感人数
lemda=0.1;%接触率
t=1:365;%时间
for i=1:(size(t,2)-1)
    %I(1+i)=I(i)+I(i)*(N-I(i))*lemda/N;
    I(1+i)=N./(1+(N/I(1) - 1).*exp(-lemda*i));
    S(1+i)=N-I(1+i);
end
plot(t,I,t,S)
xlabel('时间')
ylabel('人数')
legend('患病者','易感者')
title('SI传染病模型')